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Abstract 



We introduce a model of a nonlinear double-barrier structure, to describe 
in a simple way the effects of electron-electron scattering while remaining 
analytically tractable. The model is based on a generalized effective-mass 
equation where a nonlinear local field interaction is introduced to account for 
those inelastic scattering phenomena. Resonance peaks seen in the transmis- 
sion coefficient spectra for the linear case appear shifted to higher energies 
depending on the magnitude of the nonlinear coupling. Our results are in 
good agreement with self-consistent solutions of the Schrodinger and Poisson 
equations. The calculation procedure is seen to be very fast, which makes our 
technique a good candidate for rapid approximate analysis of these structures. 
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Resonant tunneling (RT) through double-barrier structures (DBS) make these systems 
very promising candidates for a new generation of ultra-high speed electronic devices. For 
instance, a GaAs-Gai-^Al^As DBS operating at THz frequencies has been reported in the 
literature. |I] The basic reason for RT to arise in DBS is a quantum phenomenon whose 
fundamental characteristics are by now well understood: There exists a dramatic increase 
of the electron transmittivity whenever the energy of the incident electron is close to one of 
the unoccupied quasi-bound-states inside the well. In practice, a bias voltage is applied 
to shift the energy of this quasi-bound-state of nonzero width so that its center matches the 
Fermi level. Consequently, the j — V characteristics present negative differential resistance 
(NDR). 

In actual samples, however, the situation is much more complex than this simple picture. 
This is so even in good-quality heterostructures, when scattering by dislocations or surface 
roughness is negligible. In particular, inelastic scattering is always present in real devices. An 
example of inelastic scattering events is electron-electron interaction, in which the energy of 
the tunneling electron changes and the phase memory is lost. The influence of such many- 
body effects on DBS has recently attracted considerable attention. Even with its rather 
satisfactory degree of success, many-body calculations have difficulties that, in some cases, 
may complicate the interpretation of the underlying physical processes. Recently Presilla 
et al. || pointed out the possibility of nonlocal effective nonlinearities due to many-body 
interactions in electron transport through DBS. Several results have been obtained using 
this mean field analysis such as quantum chaos [[| and nonlinear oscillations. |5||| Loosely 
speaking, this kind of treatment could be regarded as similar to Hartree-Fock and other 
self-consistent techniques, which substitute many-body interactions by a nonlinear effective 
potential. It has to be stressed, however, that the nonlocal interaction might not be the 
most suitable one in many contexts (e.g., when wells are wide) because it does not take into 
account the spatial variation of the effective potential. In any event what appears to be 
clear is that RT in nanostructures is intrinsically a (perhaps weakly) non-linear phenomena 
as have been shown in many different situations, such as the nonlinearity dependence of 
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the lock-in voltage with the number of tunneling electrons obtained with a Hartree-Fock 
approximation by Wang et al. [[?]] 

In this letter we present a model where the nonlinear interaction is driven by a local 
field instead of the mean-field approach used by Presilla and co-workers. || To describe our 
model we have chosen a DBS under an applied electric field. The thickness of the whole 
structure is L and the thickness of the well is d. The barriers are assumed to be of the 
same thickness (symmetric case) but as will be evident below this is not a restriction of our 
approach. The structure is embedded in a highly doped material acting as contact, so that 
the electric field is applied only in the DBS. We focus on electron states close to the bandgap 
and thus we can neglect nonparabolicity effects hereafter. Then the one-band effective-mass 
framework is completely justified to get accurate results. For the sake of simplicity, we will 
further assume that the electron effective-mass m* is constant through the whole structure. 
This hypothesis is related to the fact that for the time being, we are not interested in high 
quantitative accuracy, although we note that the spatial dependence of the effective mass 
can be taken into account if necessary. 

The generalized effective-mass equation is given by [we use units such that energies 
are measured in effective Rydberg (Ry*) and lengths in effective Bohr radius (a*), being 
1 R y * = 5.5 meV and 1 a* = 100 A in GaAs] 

-Mz) + \V Eff (z)-E\ ij(z)=0, (1) 

where Ve//(z) is the potential term which we discuss below. The DBS can be regarded 
as an effective medium which reacts to the presence of the tunneling electron, leading to a 
feedback mechanism by which inelastic scattering processes change the RT characteristics of 
the device. It thus follows that VEff(z) must contain nonlinear terms if it is to summarize the 
medium reaction which comes from the electron-electron and electron-lattice interactions. 
This term will be a generalized functional of the electronic probability density J r (\ijj(z)\ 2 ). 
We can now expand this functional as power series in |z/>(,2)| 2 , and neglecting higher order 
(which implies an assumption of weak nonlinearity or weak electron-electron scattering) 
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contributions we postulate that the potential in Eq. ([!]) has the form 



V Eff (z)=V X {z) + u\4>{z)\ 2 -eFz, 



(2) 



where V is the conduction band-offset at the interfaces, F is the electric field applied along 
the growth direction, and xi z ) is t ne characteristic function of the barriers, 



and all the nonlinear physics is contained in the coefficient a which we discuss below. 

There are two factors that configure the medium nonlinear response to the tunneling 
electron. First, it goes without saying that there are repulsive electron-electron Coulomb 
interactions, which should enter the effective potential with a positive nonlinearity, i.e., 
the energy is increased by local charge accumulations, leading to a positive sign for a. 
Therefore intuitively, a negative sign for a does not appear to be realistic. Nevertheless, 
we consider also that situation because in some materials, like polar semiconductors, the 
electron polarizes the surrounding medium creating a local, positive charge density. Hence 
the electron reacts to this polarization and experiences an attractive potential. This happens, 
for instance, in the polaron problem in the weak coupling limit, which becomes valid in most 
semiconductors, and where it can be seen that the lowest band energy state decreases. || 

We now work starting from Eq. (|l|) with the definition in Eq. (|2|) to cast our equations 
in a more tractable form. For simplicity, and because we are interested in intrinsic DBS 
features, we consider that the contacts in which the structure is embedded behave linearly. 
Therefore, the solution of Eq. (|I|) is a linear combination of traveling waves. As usual in 
scattering problems, we assume an electron incident from the left and define the reflection, 
r, and transmission, t, amplitudes by the relationships 



1, if < z < (L — d)/2 




(3) 



0, otherwise, 
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where k 2 , = E, k\ = E + eEL, and A is the incident wave amplitude. Now we define 
i/j(z) = (A\t\y/k~L)(j)(z), a = aki\At\ 2 . Notice that a depends on the amplitudes of the 
incoming and outgoing waves, which will be relevant later. Using Eq. ([[]) we get 

-<l> zz (z) + \y E ff(z)-E\<f>{z)=Q > (5) 

where now (Ef) read as 



V Eff (z) = V 



X (z) + a\<P(z)\ 2 -eFz. (6) 



To solve the scattering problem in the DBS we develop a similar approach to that given in 
Ref. 0. Since <f>(z) is a complex function, we write 4>(z) = q(z) exp[i r y(z)], where q(z) and 
7(2) are real functions. Inserting this factorization in Eq. (|^) we have r ) z {z) = q~ 2 (z) and 

- q zz {z) + + [VqxO) - eFz - E\ q{z) + 

q-\z) 

+V aq 3 (z) = 0. (7) 

This nonlinear differential equation must be supplemented by appropriate boundary con- 
ditions. However, using Eq. (|]) this problem can be converted into a initial conditions 
equation. In fact, it is straightforward to prove that 

q{L) = kl l/ \ q z (L)=0, (8) 

and that the transmission coefficient is given by 

4M 2 (0) 



(9) 



1 + 2M 2 (0) + A; V(0) + g 2 (0)g 2 2 (0) ' 
Hence, we can integrate numerically (0) with initial conditions (H) backwards, from z = L 
up to z = 0, to obtain g(0) and q z (0), thus computing the transmission coefficient for given 
nonlinear coupling a, incoming energy E and applied voltage V = FL. 

Once the transmission coefficient has been computed, and recalling that contacts are 
linear media, the tunneling current density at a given temperature T for the DBS can be 
calculated within the stationary-state model from 



rr)* pic D r P POO 

j{V) = -^r- j o r(E, V)N(E, V) dE, (10a) 

where N(E, V) accounts for the occupation of states to both sides of the device, according 
to the Fermi distribution function, and it is given by 

V ! \l + exp[(E F - E -eV)/k B T}J ' y 1 

where kg is the Boltzmann constant. 

In our calculations we have considered a GaAs-Gao.65Alo.35As double-barrier structure 
with L = 3d = 150 A. The conduction-band offset is Vq = 250 meV. In the absence of 
applied electric field and nonlinearities, there exists a single, very narrow resonance with 
t ~ 1 below the top of the barrier, with an energy of ~ 81 meV, and hence the well supports 
a single quasi-bound state. Figure 1 shows the transmission coefficient as a function of the 
incoming energy for different values of the nonlinear coupling a (a) and 10~ 4 , (b) 10 -3 
and (c) 10~ 2 , at zero bias. Insets show the effective potential Ve// at the energy marked 
with an arrow. It is clear that the resonances are shifted to energies higher than in the 
noninteracting case. The shift is produced by the accumulation of charge in the well as 
shown in the figures of Vb//- These results are in very good agreement with self-consistent 
calculations. [0,|n|] Note that Veh reproduce the charge accumulation in the barriers, close 
to the heterojunctions, and in the center of the well, obtained in the Hartree approximation. 



| 12|1 For completeness, we consider also the case when lattice lattice polarization effects are 
stronger than electron-electron interactions, thus leading to a negative a. This regime may 
give rise to some kind of unphysical "superconducting instability" but for small values of the 
nonlinear coefficients we have never found that problem. Figure 2 shows the transmission 
coefficient as a function of the incoming energy for different values of the nonlinear coupling 
a (a) and -10 -4 , (b) -10 -3 and (c) -10~ 2 , at zero bias. Now the resonances are shifted to 
lower energies. Also we show that for some values of the nonlinear coefficient there appears 
another quasi-bound state for lower energies. Looking at the effective potential we conclude 
that the polarization of the lattice produce a deeper well and consequently the possibility 
of new bound states, thus shifting or even doubling the original resonance. 



To apply our calculations to obtain measurable quantities, we now discuss the j — V 
characteristics, as computed from fllQp . When voltage is applied, the energy of the quasi- 
bound state level is lowered and a strong enhancement of the current arises whenever the 
Fermi level matches this resonance, thus leading to the well-known RT phenomenon. Typical 
results are shown in Fig. ^ for the particular values (a) a = 10 -3 and (b) — 10~ 3 . The 
temperature was 77 K and Fermi energy was (a) Ep = 83 meV and (b) 3 meV. In both 
cases we see clear NDR signatures. For the repulsive interaction we have a NDR peak shifted 
to higher energies (note the large Fermi energy) than in the linear case as expected in view 
of Fig |l| The attractive self-interaction causes the occurrence of a second peak at a lower 
voltage clearly related to the sideband in the transmission coefficient [Fig @(b)]. A study of 
bistability || may be carried out whitin our model; this kind of results will be obtained in 
a complete dynamical study of the problem that is in progress. 

We have demonstrated that the effective nonlinear interaction we have introduced cap- 
tures the essential physics of some inelastic effects in RT in a very simple way. The virtue 
of such an approach is that it allows to gain insight on the features of DBS without the 
burden of intensive computations providing, in an inexpensive way, a qualitative picture of 
what is to be expected in particular devices. Our model introduces a local field interaction 
that reproduce the effective potential in a more realistic way than mean-field approximation 
proposed in Ref. |Q. Our model can be easily modified to study other heterostructures in 
an electric field, like superlattices, and including different kinds of impurities and disorder. 
An important conclusion we draw from this work is that comparison of the result of current 
measurements with theoretical expectations allows to conclude which effect or group effects 
(electron-electron interaction or lattice polarization) is responsible for deviations from lin- 
earity in any specific material, by just looking at the shifting of the peaks or the appearance 
of new peaks. Also, it has to be stressed that a mere quantitative comparison between the 
present work and our idea of effective potentials to Hartree results on the same structure is 
needed, and we plan to elaborate further on this point in the new future. In addition to re- 
alize the study of the interplay between disorder and nonlinearity in this kind of models will 
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be very interesting because of the possibility of finding derealization, driven by nonlinear 



effects. 13 
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FIGURES 

FIG. 1. Transmission coefficient r as a function of the electron energy at zero bias for (a) 
a = 10 -4 , (b) 10 -3 , and (c) 10 -2 . For comparison, dashed line indicates in (a) the result for a = 0. 
Insets show the effective potential Ve// at the energy marked with an arrow. 

FIG. 2. Transmission coefficient r as a function of the electron energy at zero bias for (a) 
a = — 10~ 4 , (b) — 10~ 3 , and (c) — 10~ 2 . For comparison, dashed line indicates in (a) the result for 
a = 0. Insets show the effective potential Vg// at the energy marked with an arrow. Left inset in 
(c) shows an enlarged view of the lower resonant peak. 

FIG. 3. Computed j - V characteristics for T = 77K, with (a) a = 10~ 3 ,E F = 83 meV and 
(b) a = -1(T 3 , Ep = 3 meV. 



10 



0.8 



0.4 



0.0 
0.8 




0.4 



0.0 



0.8 



0.4 



0.0 



(b) 



250 



% 150 



> 50 



-50 



250 



% 150 - 



> 50 



-50 



-25 75 175 

z (Angstroms) 



(C) 



-25 75 175 

z (Angstroms) 




60 80 100 120 140 
Energy (meV) 



160 



(a) 



4- 



0.8 
0.6 
0.4 
0.2 
0.0 



250 



% 150 



> 50 



-50 



- 25 75 175J 
z (Angstroms) 



300 
200 
100 




•100 



(b) 



-25 75 175J 
z (Angstroms) 



300 

> 

6, 

> a -300 



-600 



(c) 



-25 75 175J 
z (Angstroms) 







40 80 120 
Energy (meV) 



160 



